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SUMMARY 


The principles of transfer function analysis have been applied to a passive 
optical heterodyne receiver to obtain the modulation transfer function (MTF). 

MTF calculations have been performed based on an optical platform which is imag- 
ing vertically varying profiles at worst-case shuttle orbit altitudes. An anal- 
ysis of the derogatory effects of sampling (aliasing) and central obscurations 
on both resolution and heterodyne efficiency is given. 


INTRODUCTION 

One measure of performance of an optical imaging system is its ability to 
reproduce an object distribution with sufficient signal-to-noise ratio and reso- 
lution so as to make the information contained within the image useful. Gener- 
ally, such a system may be characterized by its optical transfer function (OTF) 
or, in certain systems by the modulation transfer function (MTF). (See ref. 1). 

For conventional imaging systems using either coherent or incoherent illu- 
mination, one usually assumes linearity in the imaging process so that the cas- 
cading property of transfer function analysis applied (ref. 2). Under this 
assumption, the MTF's of the individual subsystems (i.e., optics, detector, elec- 
tronics, etc.) can be multiplied to give the overall transfer function. 

In this report, the principles of transfer function analysis have been 
applied to a passive optical heterodyne receiver which is assumed to be imaging 
vertically varying spatial profiles at worst-case shuttle orbit altitudes. 

Results of the analysis show some interesting departures from the properties 
described above; namely, that the cascading property must be carefully applied 
and that optical receivers having obscurations, such as Cassegrains, are not 
optimum for heterodyne-type detection. 



SYMBOLS AND ABBREVIATIONS 


Ag source geometrical shape factor 

"ig Fourier transform of *Ag, m^ 

^HIF effective mixing bandwidth, Hz 

Bsif effective shot noise bandwidth, Hz 

D;i diameter of receiver aperture, m 

D.C. direct current 

dj^ distance of image from receiver aperture, m 

Eg amplitude of geometrical image signal field, (W/m^)^/^ 

Eq amplitude of local oscillator field, (W/m^)^/^ 

Eg amplitude of signal field in detector plane, (W/m^)^/^ 

e electron charge, C 

f intermediate frequency, Hg 

Gdet Fourier transform of detector aperture function, m^ 

Gg Fourier spectrum of geometrical image spectral radiance, W/str/Hg 

Gfj normalized heterodyne transfer function 

Gq Fourier spectrum of local oscillator field modulated by detector 

pupil function, 

H coherent optical transfer function 

H(j photodetector transfer function 

Hhet total electric heterodyne transfer function 

^ifa I.F. amplifier and filter transfer function 

Hlp low-pass filter transfer function 

H^ transfer function representing carrier diffusion and transit 

time effects 

Hgq square-law detector transfer function 

H-j- tangent height 
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Htr 

h 

h(r) 

Id 

4et 

Iif(IiF^ 

Ilo 

Io»Isync 

Is 

Isync 

It 

I.F. 

I.F.O.V. 

K 

Kc 

Kob 

Kx.Ky 

K0,K4, 

KT 

LHS 

L.O. 

MTF 

N 

OTF 


photodetector shot noise transfer function 

Planck’s constant, J(sec) 

point spread function, 

photodetector current, A 

mean-square heterodyne current power, 

I.F. amplifier output current (current power); A, (A^) 

photodetector current induced by L.O. field, A 

synchronous detector output current, A 

photodetector current induced by signal field, A 

spatial Fourier component of synchronous detector current, m^A 

low-pass filter output current, A 

intermediate frequency, Hz 

geometrical instantaneous field of view 

spatial frequency vector, m~^ 

coherent cut-off frequency, cycles/rad 

coherent cut-on frequency due to central obscuration, cycles/rad 
rectangular components of K, m“^ 

angular frequency variables corresponding to K^jKy, cycles/rad 

product of Boltzman's constant and source temperature, J 

detector edge length, m 

laser heterodyne spectrometer 

local oscillator 

modulation transfer function 

shot noise current power due to L.O., normalization factor, A^ 
optical transfer function 

geometrical image spectral radiance, W/m^/sr/Hz 



Pq total L.O. power absorbed by detector, W 

R orbital height, m 

r,^ detector plane spatial vectors, m 

S/N signal-to-noise ratio in shot noise limit 

T optical transmission factor 

"^det aperture function of detector geometry 

T* conjugate of L.O, field modulated by detector aperture, (W/m2)1/2 

t time , sec 

Vq vertical component of orbital velocity, m/sec 

x,y rectangular components of m 

z distance of object from receiver aperture, m 

5(s) Dirac delta function, m“^ 

n quantum efficiency 

%et heterodyne quantum efficiency 

0,4^ angular coordinates corresponding to x,y; rad 

0p instantaneous field of view, rad 

X optical wavelength, m 

V optical signal frequency, Hz 

local oscillator frequency, Hz 

T effective integration time of low-pass filter, sec 

X efficiency factor for uniform extended source 

X, ^ efficiency factor for arbitrary source distribution 

het 

w radian frequency, rad/sec 

< > average over ensemble of signal fields 

0 convolution operation 
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Superscript: 


complex conjugate operation 


THEORETICAL ANALYSIS 
Imaging Considerations 

Consider an optical receiver which is imaging an object amplitude distribu- 
tion as shown in figure 1 . Using scalar diffraction theory, the signal ampli- 
tude Eg in the detector plane r is given by (ref. 3) 

S +oo 

d^s Eg(^)h(r - (1) 

-CO 

where h(r) is the impulse response of the imaging optics and Eg(r) is the 
amplitude of the geometrical image of the object. The shift invariance of 
h(r) can be justified for the heterodyne applications discussed here by a care- 
ful examination of the various phase factors appearing in the impulse response 
function. 

For the mixing of two deterministic optical beams in an ideal detector, the 
mean-square heterodyne current power at the difference frequency, f = |v - Vo|, 
is (ref. 4) 
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het 
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d^r E*(r) E (r) 
o — s — 
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where EQ(r) is the local oscillator amplitude distribution in the detector 
plane, ri is the quantum efficiency, e is the electronic charge, and hv 
is the photon energy. A simple classical approach is taken to obtain the 
correct expression from which the spatial frequency analysis may begin. We 
recognize that the geometrical image field Eg(s,v) is a stochastic process 
which we synthesize by discrete frequency components with random phases. Now, 
for a deterministic L.O. field and a quasi-monochromatic optical signal, equa- 
tion (1) and the generalization of equation (2) combine to give 
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where < > represents an average over the ensemble of signal fields. It is 
assumed that the source of the image field on the detector (i.e., the sun) is 
spatially incoherent. The appropriate substitutions are 



I 


2 

het 




where Pg represents the image spectral radiance at the detector plane in 

W/m^/sr/Hz and .(f) is the current spectral power density in A^/Hz. 

het 

Equation (3) becomes 


f) ]|t*(s) ® h(-s) 1^ 

where f ^ 0 and we note that the L.O. mixes with the signal field components 
at Vq + f and Vq - f. Here, we have expressed the detector overlap integral 
of the L.O. field and the impulse response function as 

+00 

1 Tj . (r)E*(r)h(r - s)d^r = T*(s) Q h(-s) 

J det — o — — — o— — 

— oo 

where is the aperture function of the detector geometry and the pro- 

duct T(jet is simply that portion of the L.O. that is trans- 
mitted by the detector aperture. 

Referring to the detector scheme of figure 2, the output current from the 

J + Pg (s , v^-f) ] | T* (s) 

(f>0) 


synchronous detector is 
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where is the total heterodyne transfer function defined by 






and T is the optical transmission factor. The various contributions to 
Hjiet are: the signal/L.O. mixing transfer function Hjjj, representing carrier 

diffusion and transit time effects in the detector; the photodetector transfer 
function comprised of contributions due to capacitance, resistance, and 

inductance; and the I.F. (intermediate frequency) amplifier and filter transfer 
function The square-law detector is assumed to have a unity transfer 

function (Hgq = 1). The shot-noise transfer function is due only to 

transit time effects as opposed to Hjjj. 


Heterodyne Transfer Function 

In this paper, we are interested in the spatial frequency response of the 
heterodyne receiver to a vertically varying object profile as shown in figure 1 . 
Other than a knowledge of the total I.F. power, this is different from the I.F. 
considerations discussed previously. To obtain the spatial frequency response 
and, ultimately, the system modulation transfer function (MTF), we assume the 
object scene radiance (and, consequently, the image scene) is linearly trans- 
lated owing to motion of the optical receiver (e.g., an orbiting platform). 

This induces a translation of the image coordinates by an amount 


s ->■ s - r(t) 


and 


Pg(S, + f) ^ Pg(s 




Further, we define 

Pg(-r) = J df|H^^^(f)|2[p^(r, + f) + Pg(r, - f ) ] 

(f>0) 

Since the impulse response h(r) is invariant, we have the output current from 
the synchronous detector as 


I 

sync 


(t) 


I [r(t) ] = 
sync — 


2Te^n^X^ 

(hv)^ 
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(4) 
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Equation (4) is of the form of a convolution 


P_(r) ® |T*(r) Sh(-r )|2 
sync g o 


by the 


Decomposition of Isync ^ 
e Fourier transformation 


into its spatial frequency components is obtained 


I (K) 
sync — 


oo 

= I 


-a 27 rK*r ^ , x 

e I (r) d r 

sync — 


By use of the convolution theorem, we have 


1 (K) = ^ G (K) [g*(-K)H(-K) a G (K)H*(K)] (5) 

~ (hv) ^ g-o - - o- - 


where K is the spatial frequency vector variable defined by its rectangular 
components (K^jKy), Gg(K) is the object, or more specifically, the geometrical 
image spectrum, i® the detector pupil function modulated L.O. spectrum, 

and H(K) is the coherent transfer function of the system (ref. 3). 

Equation (5) illustrates the departure of the transfer function obtained 
in a heterodyne system with that obtained in conventional imaging systems. 
Remembering that the coherent transfer function H(K) is equal to the pupil 
function of the optical receiver (with a suitable change in variables) (ref. 3), 
the conventional optical transfer function is proportional to 

Gdet<K) [h(-k) SH*(K)] 


where G^jet(K) is the Fourier transform of the detector aperture function, 
Tdet^^)* equation (5), we see, however, that H(K) is modified by the 
spectrum of the L. 0/detector combination Gq(K). The normalized convolution 
of the product Go(K)H*(K) with its negative argument complex conjugate is 
defined as the heterodyne transfer function Gjj. Functionally, then, we define 
a normalized heterodyne transfer function by 



or 


°H® - 
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d\'G (K')H*(K')G*(K’-K)H(K'-K) 
o — — o — — — 


d^K' |g (K’) I^Ih(K') 


Heterodyne Efficiency Factor 


The denominator of equation (6) indicates that the product Gq(K)H(K) 
represents the optics/L.O. detector amplitude spectrum that is transferred to 
the detector. Using Parseval's theorem, the integral 


5d2K| Go(K) I 2| H(K) I 2 is thus the power available 
total L.O. -detector power of Pg H^d^K] Gg(K) | ^. 


for heterodyning out of a 
For a uniform extended source, 


we may thus define an efficiency factor 


Jd^KlG^(K) |^|H(_K) 

X = 1 1 (7) 

Jd\|G^(K) 1^ 


With this definition, the current spectrum for the synchronous detector (equa- 
tion (5)) becomes 


Wnc® -^^^V<=H®\® 


( 8 ) 


Equation (8) may be related to a more conventional form of heterodyne 
efficiency found in the literature (ref. 5). The synchronous detector current 
is the inverse Fourier transform of equation (8); that is, 
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For a stationary scene (i.e., before translation of the image coordinates), we 
have r = 0 so that 
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sync 


( 0 ) 
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—09 


i (K) d^K 
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Substitution of equation (8) into equation (9) and assuming a blackbody source 
of geometrical shape factor Ag(r) and radiance Pg(r), we have 


PgW - 




\(£> 


and 


I (r = 0) 
sync — 




+QO 

J d^K Ag(K)G^(K) 

— OO 


( 10 ) 


where mixing occurs over an effective bandwidth 2Bfjjp centered at the L.O. 
frequency, a polarization loss factor of 0.5 is included, and Ag(K) is the 
Fourier transform of Ag(r). The integral portion of equation (TO) has the 
form of a throughput (i.e., that portion of the image passed by the heterodyne 
transfer function) . The product of X and this integral is an efficiency 


^het 
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yj 


d^K 


Ag(K)Gjj(K) 
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SO that equation (10) becomes 
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we define a shot-noise level due to the L.O. by 


If now, 
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2ne^ 

(hv) ^o^SIF 


then the signal- to-noise ratio in the shot-noise limit becomes 
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where we have defined 
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het 


>^het ^ 
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as the heterodyne quantum efficiency and Bsip is the effective shot noise 
bandwidth. Note that for a uniform extended source, Ag(K) = 6(K) and equa- 
tion (11) reduces to Xhet = X and Hhet " BX* In this case, the efficiency 
factor x» which we have defined in equation (7), is equivalent (to within the 
D.C. quantum efficiency q) to the heterodyne quantum efficiency, hhet* f*onnd 
in the literature (ref. 5). 


System Transfer Function 

Results from the previous section may now be used to calculate the system 
transfer function, including the low-pass filter (see fig. 2), for the specific 
example of imaging a one-dimensional object through an optical receiver which 
has rectangular symmetry. This has some physical significance since the reso- 
lution elements of interest in an orbiting heterodyne receiver are vertically 
varying stratospheric layers. To avoid scaling difficulties in the calcula- 
tions, we will use angular coordinates defined by (see fig. 1) 



( 12 ) 

K0 = Kydi K(j) = Kydi (cycles/rad) 
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and 


0F = ;s- (rad) 


where 0p is the geometrical instantaneous field of view (I.F.O.V.) of the 
optical receiver. 

Using equations (12) and the one-dimensional geometry, equation (6) becomes 




Kq k 

[sincOpK )RECT(— ) ] S [sine ( ) RECT(— ) ] 
c 

{[sincOpKglRECTC^) ] 0 [sinc(0pKg)BECT(— ) ]}^ 
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(13) 


where 


_ sin i:X , ,X , _ 

sxnc(X) = — and RECT(— ) = 


1 , X < t 


0 , lx| > t 


Equation (13) assumes a plane-wave local oscillator incident of the detec- 
tor so that the detector/L.O. transfer function becomes simply the Fourier 
transform of the detector aperture. Further, the coherent transfer function 
for the optics is the pupil function (rectangular in shape) having a coherent 
cutoff frequency of = Dj^/2X, where is the diameter of the receiver 

aperture and X the wavelength. This convolution process is shown in figure 3. 

Equation (13) along with equation (8) gives the system transfer function 
up to the low-pass filter. Expressing G{j(K0;Kq,0f) and X(Kc»9f^ (7)) 

in integral form, we have 
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( 14 ) 



and 


+00 


x(Ke.0p) 
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9 K 

sine (e^Kg) SECT (— ) dKg 


( 15 ) 


Finally, inclusion of the low-pass filter transfer function H^p (i.e., 
the integration time, t) which is modeled as a running mean integrator, we have 




where Vq is the vertical component of the orbital velocity, z the receiver- 
object distance, and x the integration time. The total transfer function is 
then the modulation transfer function 


Equations (14) and (15) can be evaluated in terms of tabulated functions 
yielding the following relations which will be used for computational purposes 
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In the above equations, the functions 
(ref. 6) 

. 

s. (X) = \ j " dt and C. (X) 
1 ] t xn 

•^o 


RESULTS AND DISCUSSION 
MTF Calculations 

Equations (T6) and (17) may now be evaluated for some specific parameter 
values which are applicable to the optical receiver in a spacelab type of 
scenario. A worst-case set of orbital values would be obtained for the receiver 
platform at an orbital height of R □ ilOO Km and a tangent height Ht of 10 km. 
At these values, we assume that the receiver is operating in a solar occulta- 
tion mode where the sunrise or sunset velocity due to orbital motion is Vq = 

2 km/ sec. The MTF and X calculations (equations (16) and (17)) will be done 
for an I.F.O.V. of 0p = 0.5 mrad, an equivalent optical receiver aperture of 
D;^ o 50.8 mm, and values of integration time of t = 0.2 sec and 0.4 sec. Fur- 
ther, the value of D^ = 50.8 mm at a wavelength of X z 11.1 52 pm (HNO3 line) 
corresponds to an optics cutoff frequency of Kq = D^/ 2X z 2278 cycles/rad. 

These parameters are compatible with the values for an LHS type experiment using 
a tunable diode laser as the L.O. and associated optics for coupling this type 
radiation to a detector having the required time- frequency response (ref. 7). 

The results of the calculations are shown in figure 4. It can be shown 
that, for values of t greater than roughly 0.2 sec., the optical (heterodyne) 
transfer function dominates the MTF; and, for t somewhat less than 0.4 sec., 
the low-pass filter is the dominant frequency limiting factor. Note that the 
angular frequency values can be converted to linear spatial frequency (cycles/ 
km) by the relationships of equation (12) by appropriately scaling image and 
object space by the ratio of image distance di to object distance z. For 
the orbital values assumed, z = 2262 km and consequently a value of 2262 
cycles/rad corresponds to an object spatial frequency of 1 cycle/km. Examina- 
tion of the MTF curves shows that resolutions of the order of 1.5 to 2.0 km may 
be expected for the various integration times. 


Sj[(X) are defined as 



Efficiency Calculations 

The efficiency factor (heterodyne efficiency) given by equation (17) is 
shown plotted in figure 5 for the example of the plane-wave L.O.. Two geome- 
tries are shown: rectangular optics (as has been previously assumed) and cir- 

cular optics adjusted for equal optics and detector areas. The significance 
of the abscissa (2 tt 6 pKq) relative to heterodyne efficiency becomes apparent when 


14 



it is noted that at the value of 2TT0pKQ 7.7 corresponds approximately to 
an image (sun) size filling the detector of one Airy disk of the receiver 
aperture. In this region, the efficiency is in excess of 80 percent. 


Sampling Error 

The calculations shown plotted in figure 5 do not include any sampling 
errors which may occur. Suppose we sample the output of the low-pass filter, 
which has been modeled as a running mean integrator, at a rate of the inverse 
of the integration time. This is equivalent to a sampled mean integration 
scheme. Under this constraint, it may be shown that for certain values of T 
the signal is undersampled. This results in an aliasing or foldover error 
which can be significant relative to the desired signal. For example, shown 
in figure 6 is the total MTF for the values of orbital and system parameters 
previously stated. Two integration times are considered: x = 0.2 sec and 

X = 0.4 sec. If we define the sampling error as the ratio of the ” foldover” 
amplitude on the MTF plot to the amplitude of the MTF itself (i.e., a white 
signal spectrum), we see that the error for 0.4 sec. is approximately 40 per- 
cent at 0.5 cycle/km frequency and considerably worse for higher values of 
spatial frequency. Conversely, for x = 0.2 sec and the correspondingly 
higher sampling rate, the sampling error is negligible. 


Heterodyne Receivers With Obscurations 

Telescopes having central obscurations, such as Cassegrains, are often 
used for imaging a source. If this type receiver is used as a collector for 
heterodyne-type detection, one needs to compare the efficiency x and the 
heterodyne transfer function with that obtained for no obscuration case. 

In figure 7, we consider the effects of receiver apertures having obscura- 
tion ratios of 0 and 20 percent for Op = 0.2 and 0.5 mrad. Note the enhance- 
ment of response in the 2000 cycles/rad region at the expense of that near 1000 
cycles/rad for 20 percent obscuration and Op = 0.5 mrad. The effects of 
obscurations are more pronounced for square as opposed to circular geometries. 

An unobscured conventional MTF discussed earlier is plotted for Op = 0.5 mrad 
showing a somewhat reduced frequency response characteristic from the heterodyne 
MTF. For a smaller detector (Op = 0.2 mrad), a 20 percent obscuration tends 
to assume the shape of a conventional MTF. 

In figure 5, we assume a receiver aperture having obscuration ratios of 
20 and 50 percent. For the values of Op and used earlier, 27r0pKQ = 7,2, 

and comparison of the various efficiency curves at this value shows striking 
differences. For 50 percent obscuration case, one sees that the heterodyne 
efficiency is virtually zero while for 20 percent obscuration a relative effi- 
ciency of slightly greater than 0.2 is achieved. This compares with a value 
of greater than 0.8 when unobscured. Note further, there is a "peaking” of the 
efficiency curves for obscured systems. The rule-of-thumb requirement of one 
Airy disk on the detector for "good" efficiency no longer holds, but rather the 
source image needs to be less than this value to achieve the maximum efficiency 
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for that particular system. The effect can be explained on the basis of the 
overlap integral (eq. (2)) of the L.O. field and signal field distributions. 
Thus, the diffracted field due to the central obscuration is out of phase with 
that of the primary diffracted field and, as the size of the detector and/or 
optics increases, the cancellation tends to be more complete. 


CONCLUDING REMARKS 

The analysis of a passive heterodyne receiver with respect to its imaging 
performance (transfer function) and its heterodyne efficiency shows some 
interesting departures from the results which are obtained in strictly coherent 
or incoherent imaging systems. For example, the cascading property of MTF 
analysis must be carefully applied, since the coherent transfer function of the 
optical receiver and that due to the L.O. -detector combination are not separable 
but are related by the convolution of their products. Application of these 
results to a spacelab-type optical heterodyne receiver (LHS) shows that resolu- 
tions of the order of 1.5 to 2.0 km are possible for worst-case type orbital 
scenarios . 

Further, comparison of obscured-type receivers (e.g., Cassegrains) with 
unobscured receivers shows that both resolution and efficiency are severely 
degraded in an obscured-type receiver and, consequently, should not be used 
for a passive heterodyne detection scheme. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 24 , 1980 
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Figure 2.- System transfer functions. 
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Figure 5 



- Efficiency factor, X/ versus system parameter 27r0pK^ for various 

receiver geometries. 



Figure 6.~ Aliasing error for 2.5-Hz sampling rate at worst-case shuttle 

orbit. 



TRANSFER FUNCTION 



Figure 7.~ Comparison of transfer functions for various geometries. Dashed 
curve represents conventional MTF with 0.5 mrad I.F.O.V. 
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lutions of the order of 1.5 to 2.0 km are possible for worst-case type orbital 
scenarios. In addition to an analysis of the derogatory effects of sampling 
(aliasing), a comparison of obscured-type receivers (e.g., Cassegrains) with 
unobscured receivers shows that both resolution and efficiency are severely 
degraded in an obscured-type receiver and, consequently, should not be used for 
passive heterodyne detection schemes. 
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